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Abstract-In this paper, algorithms to estimate narrowband to 
broadband albedo (NBBA) of snow using Advanced Wide Field 
Sensor (AWiFS) and Moderate Resolution Imaging Spectro- 
radiometer (MODIS) sensor images have been developed for 
Indian western Himalayan region. The in-situ measurements of 
spectral reflectance and transmitted spectral solar irradiance 
on snow surface by spectroradiometer have been used to 
calibrate and validate the algorithms developed. The AWiFS 
and MODIS derived snow broadband albedo via the developed 
algorithms have been validated with in-situ observations at 
three field locations. An RMSE better than 0.03 and an R 2 
value of 0.94 and 0.88 between modeled and observed albedo 
values have been obtained for AWiFS and MODIS respectively. 
The algorithm developed for MODIS images has also been 
compared with an existing global albedo model for MODIS, 
and the results from the proposed model have been found to be 
in agreement with the published model. 
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I. INTRODUCTION 

Broadband snow albedo in solar wavelength region (0.3 
|um - 3.0 |um) has been an important parameter for various 
climatic studies. It is a key input to local or basin scale 
snowmelt runoff models and is also a significant parameter 
to understand and characterize the energy fluxes at the earth 
surface It also happens to be a pointer to understand 
many physical aspects of the snow surface, such as snow 
depth, grain size and water content [2] . Thus, the spatial 
variability of snow surface albedo is of immense importance 
in many applications such as ecosystem dynamics, 
hydrology, geomorphology and microclimate of high relief 
environments [3] . Variability in snow albedo exerts 
significant impact on snow cover properties during the 
course of accumulation and ablation. 

Large stretch of mountainous regions in Himalaya is 
covered by permanent snow cover, which increases 
manifolds during winter season. Seasonal snow cover 
generally builds up and ablates from November to June in 
different ranges of lower and middle Himalaya [4] . Recent 
analyses of long term observations of snow precipitation 
show a declining trend of snowfall in Beas basin of 
Himalayan region [5] . The snow cover in mountainous 
region results in a cooler climate there by affecting the local 
energy balance through reduced absorption of the solar 
radiation. The study of snowmelt of the seasonal snow cover 
therefore requires an understanding of spatial and temporal 
variability of albedo in the Himalayan region. 



Satellite remote sensing offers a means to quantitatively 
estimate the surface reflectance in discrete spectral bands in 
large spatial area of inaccessible Himalayan terrain. 
However, this spectral reflectance may not be directly used 
for energy balance estimation as it requires a short wave 
integrated albedo in the reflective solar wavelength region. 
Therefore, an appropriate methodology to estimate 
broadband integrated albedo over a snow-covered area from 
spectral reflectance of snow in discrete satellite bands is 
required. The narrowband reflectance from medium (e.g. 
IRS AWiFS) and moderate (e.g. Terra MODIS) spatial 
resolution satellite data may have the potential to provide an 
estimation of integrated snow albedo using the NBBA 
algorithms at frequent time intervals due to their high 
temporal resolution. 

Reference [6] has reported an extensive work on 
conversion of narrowband to broadband albedo (NBBA) for 
a number of satellite sensors including Terra Advanced 
Space borne Thermal Emission and Reflection Radiometer 
(ASTER), NOAA Advanced Very High Resolution 
Radiometer (AVHRR), Geostationary Operational 
Environmental Satellite (GOES Imager), LANDSAT 7 
Enhanced Thematic Mapper Plus (ETM + ), Terra Multiangle 
Imaging Spectroradiometer (MISR), Polarization and 
Directionality of Earth's Reflectances (POLDER), SPOT 
VEGETATION and Terra MODIS. The listed algorithms 
have been testified under different atmospheric conditions 
and solar geometry through MODTRAN, discrete-ordinate 
radiative transfer (DISORT) and other atmospheric 
correction models. Previous studies have established 
accurate estimation of broadband albedo from the Advanced 
Very High Resolution Radiometer (AVHRR) satellite data 
in polar region [7 ' 8] . Due to improved spectral/spatial 
characteristics and measurement precision of MODIS 
images as compared to AVHRR, it is anticipated that 
accuracy of surface albedo estimation may further increase. 
Reference [9] has applied MODIS Level IB data to derive 
the albedo over Greenland. The satellite-derived surface 
albedo has been found to be within about 6% of that 
measured at the stations. Reference [1] has developed a 
prototype snow albedo algorithm that also relies on MODIS 
atmospherically corrected surface reflectance. However, it 
accounts for reflectance anisotropy using DISORT model. 
Reference [10] has developed a direct retrieval algorithm 
that links top-of-atmosphere (TOA) narrowband albedo to 
land surface broadband albedo using a feed-forward neural 
network. The algorithm, although, has been designed for 
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estimating albedo over any land cover types but appears to 
be particularly well designed for non-vegetated land cover. 
Moreover, it does not take into account atmospheric 
correction and other processes. Reference [11] has given a 
series of formulae for conversion of NBBA based on 
extensive radiative transfer simulations for data from 
different satellite sensors. One of the algorithms given in 
[10] has later been revised and applied to MODIS data for 
global albedo estimation [12] . 

The AWiFS data have very different spatial and spectral 
resolutions as compared to data from other remote sensing 
sensors. Therefore, the broadband albedo model developed 
for other sensors may not be directly applicable to AWiFS 
data. In the present study, an algorithm for narrowband to 
broadband albedo of snow using AWiFS data has been 
developed. The same methodology has also been extended 
to develop the algorithm for NBBA of snow using MODIS 
data. The efficacy of the algorithm for NBBA has been 
tested with the published algorithm for estimation of global 
albedo of snow using MODIS. The proposed algorithm uses 
narrow band reflectance of snow in visible (VIS), near- 
infrared (NIR) and short wave infrared (SWIR) bands. 

II. STUDY AREA DESCRIPTION 

The first study area lies between latitude 32° 12'N to 32° 
60'N and longitude 77 °E to 77°30'E and has been marked as 
ABCD on a map of India (Fig. la). 
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Fig. la) First study area marked with rectangle (ABCD) on India map 



The AWiFS image corresponding to this area is shown 
in Fig. lb. 
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Fig. lb) First study area ABCD and three field locations on AWiFS image 
for in-situ measurements of albedo for validation with satellite 



The study area includes three scientific field 
observatories, namely, Solang (latitude 32°19'27" N, 
longitude 77°09'27"E), Dhundi (latitude 32°21'22 ?, N, 
longitude 77°07'42"E) and Patsio (latitude 32°45'11"N, 
longitude 77°15'42"E) of Snow and Avalanche Study 
Establishment (SASE) in India. These observatories are 
responsible for collection of a variety of in-situ data 
including albedo measurement of snow at different time 
intervals. Solang and Dhundi observatories are located in 
the Beas basin in the Pir-Panjal Mountain range at altitudes 
of 2480 m a.s.l and 3050 m a.s.l respectively. The area is 
surrounded by forest and tree line which exists up to 3100 m 
a.s.l.. The Patsio observatory lies in the Greater Himalayan 
range at an altitude of 3880 m a.s.l and is devoid of forest. 
The two mountain ranges have different climatic and 
geomorphologic conditions. The monthly averaged 
minimum temperature during winter has been observed to 
be from -4.1°C to 2.5°C in Pir-Panjal range and from -6.4°C 
to -15.4°C in Greater Himalayan range [4] . Pir-Panjal range 
receives higher seasonal snowfall (with season average of 
10-12 m) as compared to Greater Himalayan range (with 
season average of 5-7 m). The slope in the study area varies 
from 1 to 86 degree with a mean value of 28 degree. Most 
of the slopes in this study area are southern. 

The second study area, marked as EFGH in Fig. la, lies 
between latitude 30°0'0 M N to latitude 36°0'0"N and 
longitude 73°0'0 M E to 80°0'0 M E. The MODIS image 
corresponding to this area is shown in Fig. 2. 



Latitude - 35 a 53 J 15.19" N 
Longiucte - 73° W 53-35" E 




Latitude - Stf 5 12 h 5E.SB U N 
L(mgitud*-79 D 45'<ia25- , 'E 



Fig. 2 Different Himalayan ranges corresponding to EFGH (Fig. la) on 
MODIS image used for proposed broadband albedo model comparison 
with reported global (Liang model) model 

This area, having a larger areal extent than the first study 
area, has been used for validation of the proposed algorithm 
and its evaluation vis-a-vis published global albedo model 
for MODIS. This has been done to assess the efficacy of the 
proposed algorithm for a large snow covered area having 
varied topographic and meteorological conditions. The area 
mainly encompasses three different Himalayan zones, e.g., 
Lower Himalayan or Subtropical Zone, Middle Himalayan 
or Mid Latitudinal Zone and Upper Himalayan or High 
Latitudinal Zone, as described in detail in [13]. 

III. IN-SITU MEASUREMENTS OF SPECTRAL REFLECTANCE 

Extensive in-situ measurements of spectral reflectance 
and transmitted solar irradiance on snow in the wavelength 
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range from 350 to 2500 nm have been recorded 
simultaneously at the time of satellite passes near Solang 
and Dhundi observatories. Field Spec Pro spectroradiometer 
(M/S Analytical Spectral devices, USA), shown in Fig. 3a 
has been used to collect measurements of snow at 25° field 
of view (FOV) under clear sky conditions on nearly flat 
surface. Simultaneously, solar transmitted spectral 
irradiances have also been measured using a full sky remote 
cosine receptor having a 180° FOV. These in-situ 
measurements have been used for NBBA model 
development and validation of the data acquired from 
AWiFS and MODIS at Solang and Dhundi sites. The results 
from NBBA model for AWiFS and MODIS data have also 
been validated at Patsio, where in-situ measurements of 
albedo have been collected using Albedometer at the time of 
satellite passes. Albedometer, mounted on an Automatic 
Weather Station (AWS), consists of two CMP3 
pyranometers (Kipp and Zonen) with one looking up 
towards sky and another looking down towards snow 
surface (Fig. 3b). 




Fig. 3a) Optical Spectro-radiometer and its different components at snow 
covered Dhundi location 




Fig. 3b) Albedometer (two pyranometers CMP3) for broadband albedo 
measurements mounted on Automatic Weather Station at Patsio location 



IV. SATELLITE SENSOR DATA 

Almost clear sky multi-spectral medium spatial 
resolution (56 m) and 10-bit radiometric resolution data 
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obtained from AWiFS on-board RESOURCESAT-I and 
moderate resolution (250 m and 500 m) having radiometric 
resolution 12 bit obtained from MODIS on-board Terra 
have been used in this study. The salient characteristics of 
AWiFS and MODIS sensors have been given in the Table I 
and Tablell respectively. Twelve AWiFS images and nine 
MODIS images (Table III) of different dates have been used 
for NBBA model development and its validation using in- 
situ observations. 



TABLE I SALIENT SPECIFICATIONS OF AWIFS SENSOR 



Spectral 
Bands 


Spectral 
Wavelength 
(nm) 


Maximum 
Radiance 
(m W/cm 2 /sr/um) 


Solar 
Exoatmospheric 
Spectral 
Irradiance 
(mW/cnr/um) 


B2 


520-590 


52.34 


185.3281 


B3 


620-680 


40.75 


158.042 


B4 


770-860 


28.425 


108.357 


B5 


1550-1700 


4.645 


23.786 



TABLE II SALIENT SPECIFICATIONS OF MODIS SENSOR 



Spec- 
tral 
Bands 


Spectral 
wave- 
length 

(nm) 


Spatial 
Resolu- 
tion (m) 


Radiance_ 
Scale 
(mW/cm 2 /sr/u 
m) 


Solar 

Exoatmosphe- 
ric Spectral 
Irradiance 

(mW/cnr/um) 


Bl 


620-670 


250 


0.0026144 


160.327 


B2 


841-876 


250 


0.0009926 


98.70 


B3 


459-479 


500 


0.0027612 


209.071 


B4 


545-565 


500 


0.0021087 


186.4 


B5 


1230- 
1250 


500 


0.0005568 


47.6 


B6 


1628- 
1652 


500 


0.0002572 


23.8 


B7 


2105- 
2155 


500 


0.0000787 


8.7 



TABLE III DATES OF ACQUISITION OF AWIFS AND MODIS IMAGES USED 



S.No. 


AWiFS Images 


MODIS Images 


1 


19 January 2005 


26 December 2008 


2 


21 February 2005 


31 December 2008 


3 


26 February 2005 


20 January 2009 


4 


13 March 2005 


28 January 2009 


5 


7 March 2007 


30 January 2009 


6 


8 January 2009 


26 February 2009 


7 


24 December 2008 


16 March 2009 


8 


31 January 2009 


13 April 2009 


9 


15 February 2009 


17 April 2009 


10 


25 February 2009 




11 


11 March 2009 




12 


16March 2009 





V. TOPOGRAPHY OF STUDY AREA 

A digital elevation model (DEM) of the study area 
ABCD (Fig. la) has been generated from on-screen 
digitization of contours on topographical map at 1:50000 
scale at 40 m contour interval and is shown in the Fig. 4a. A 
non-linear interpolation function has been used for DEM 
generation. The DEM has been re-sampled from 6 m to 56 
m to match the spatial resolution of AWiFS image. From 
this DEM, other topographic parameters, namely, slope (Fig. 
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4b), aspect (Fig. 4c) and local illumination image have been 
derived for terrain corrections using slope match method. 



I (b) 




Fig. 4a) Digital elevation model; b) Slope and c) Aspect of study area 
ABCD at 6m spatial resolution 

VI. METHODOLOGY 

Snow covered surfaces are usually characterized by very 
high spectral reflectance in the visible wavelength region, 
which rapidly decreases in NIR and SWIR regions. Thus, 
reflectance in narrow bands in VIS and SWIR regions of 
AWiFS and MODIS sensors appear very dominant in 
characterizing snow surface conditions. Therefore, data 
from three bands of AWiFS (Bands 2, 4 and 5) and seven 
bands of MODIS (Bandl-Band7) have been used to develop 
the algorithm for NBBA of snow. The flow chart of detailed 
methodology for NBBA is given in the Fig. 5. 
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BBA using AWiFS 
and MODIS images 



Field measurements - 

spectral reflectance 
and spectral irradiance 

V 
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V 
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Validation of BBA for AWiFS and 
MODIS with field measurements 
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Fig. 5 Flow chart showing methodology for model development of 
broadband albedo and validation 



First, image based atmospheric corrections for path 
radiance without considering topographic influences and 
Rayleigh scattering was performed for initial estimates of 
reflectance at each pixel of AWiFS and MODIS images [14, 
15] . The AWiFS reflectance image was further refined by 
applying topographic correction using slope matching 
method, which is reported to be the most suitable method 
for Himalayan region [16] . The topographic correction on 
MODIS reflectance image has not been applied due to its 
moderate spatial resolution [17] . 

All the AWiFS and MODIS images have been 
geometrically and atmospherically corrected before the 
estimation of spectral reflectance at each pixel. 

A. Geo-Referencing 

Available geocoded IRS -LIS S -III (Indian Remote 
Sensing, Linear Imaging Self Scanning) image was used as 
the master image for geo-referencing of AWiFS images. 
The MODIS images were then georeferenced to AWiFS 
images. The geo-referencing of images has been performed 
using second order polynomial transformation at root mean 
square error of less than a pixel estimated using 18-20 
ground control points (GCPs). 

B. Atmospheric Correction 

First, the digital numbers of pixels in AWiFS images 
have been converted to at- sensor spectral radiance, 
L S atA (mWcm^Sr'Vrn _1 ) using the following relationship, 



dAWiFS) 



DN, 



sat A 



+ L n 



(1) 



where DN max x and DN sat ^ are maximum and at-satellite 
DN values respectively in different AWiFS spectral bands. 
L max x and L min x are the maximum and minimum values of 
spectral radiance in different bands, as obtained from 
sensors' header information. h minX - and DN max ^= 1023. 

At satellite spectral radiance, L satA , for each pixel of 
MODIS images in different bands is calculated using the 
(2) : 

L sat/L (MODIS) = 
radiance _s c ale Sx x (DN — radiance_offset^) 

(2) 

where radiance _scales % and radiance_offset^ (zero for 
all bands) are the scale and offset parameters, respectively 
as given in Table II. 

The estimated values of at-satellite radiance at pixel 
were corrected using combination of dark object subtraction 
(DOS1 and DOS3) model [15] for atmospheric correction. It 
is a simple and efficient model for operational use in the 
absence of in-situ data of various atmospheric parameters 
such as aerosol, water vapor content, etc., at the time of 
image acquisition [18 ' 19] . The DOS1 model assumes no 
atmospheric transmittance loss and no diffused downward 
radiation at the surface [14] . DOS 3 computes up welling and 
downwelling radiation transmittances of the atmosphere to 
correct the image for path radiance resulting from the 
interaction of the electromagnetic radiation with molecules 
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and aerosols. It assumes Rayleigh scattering only. The 
optical thickness for Rayleigh scattering is obtained from 
the equation proposed by [20]. 

The values of transmittance using Rayleigh scattering 
are estimated for AWiFS and MODIS spectral bands and 
shown in the Table IV. 



TABLE IV UPWELLING AND DOWNWELLING TRANSMISSIVITY DUE TO 
RAYLEIGH SCATTERING FOR AWIFS AND MODIS SPECTRAL BANDS 



Transmiss- 
ivity 


Ban- 
dl 


Ba 
nd 

2 


Ba 
nd 

3 


Ban 
d4 


Ban 
d5 


Ban 
(16 


Ban 
d7 


AWiFS-"t v " 




0.9 
35 


0.9 
65 


0.98 
6 


0.99 
9 






AWiFS-"t z " 




0.9 
11 


0.9 
52 


0.98 



0.99 
8 






MODIS- 

"t v " 


0.965 


0.9 
89 


0.8 
77 


0.93 
6 


0.99 
7 


0.99 
9 


1.0 


MODIS-"t z " 


0.941 


0.9 
81 


0.8 
01 


0.89 
5 


0.99 
6 


0.99 
9 


1.0 



The dark objects (e.g., shadow, deep clean water) are not 
absolutely dark [21] due to atmospheric scattering effects. 
Assuming 1% surface reflectance for the dark objects, the 
path radiance was estimated using the proposed model [14 ' 22 ' 
23] . The spectral reflectance corrected for path radiance, 
under lambertian assumption, of each pixel in AWiFS and 
MODIS images is then computed using the following (3) [15 ' 

24]. 

R x = " (W " Lp)d2 (3) 

t v (E cos Z t z +E d ) 

where t v and t z are the transmittances of the atmosphere in 
the viewing and illumination directions respectively, as 
determined from the equations proposed [15 ' 20] . E is the 
exo-atmospheric spectral irradiances as given in Tables I 
and II. 6 Z is the solar zenith angle at each pixel and 
estimated using [25] d is the earth- sun distance in 
astronomical units obtained from the approach given by [26]. 
E d is the downwelling diffused radiation and has been 
assumed as zero according to [22]. L p is the path radiance 
obtained from the approach given in [22] 

C. Topographic Corrections Using Slope Match Method 

A number of image based topographic correction 
methods (e.g., Cosine correction, C-correction, Minneart 
correction, statistical etc.) have been proposed [27 ' 281 . An 
extensive study on the feasibility of different topographic 
correction methods for the Himalayan region has been 
reported [16] and recommended slope matching method more 
appropriate for the region. Slope matching method is given 
in (4) and used for topographic correction in the present 
study: 

R l = R A + (R^max - R Amin ) ( <CQS . >s ) Q ( 4 ) 

where R\ is topographically corrected spectral 
reflectance, is obtained from (3), R^ max and R^nin are 
the maximum and minimum spectral reflectance values as 
obtained from R x image, < cost > s is average value of 
illumination on south aspect, cost is illumination obtained 



from the equation proposed by [27] . C x is the normalization 
coefficient for different bands and is estimated using 
equation given by [29]. 

D. Narrow Band Reflectance (NBR) to Broadband Albedo 
(BBA ) From AWiFS and MODIS 

Reference [30] has proposed an approach to convert 
NBR to BBA from Advance Very High Resolution 
Radiometer (AVHRR) image acquired in two bands. The 
same approach has been adopted here to derive a new 
algorithm considering three spectral bands (Band 2, Band 4 
and Band 5) of AWiFS image (Table I) and seven spectral 
bands (Band 1 - Band 7) of MODIS image (Table II). 

The surface integrated albedo, a, in any wavelength 
interval, AA, defined as the ratio of the integrated reflected 
energy , £ E sAX P AX , to the integrated downward irradiance, 
Y, E sA a , at the solar zenith angle, 6 Z and can be expressed as 

a _ Iaa^saa p A a ^ 

T.AA E sAA 

where E sAA is the incoming solar energy and p AA is the 
surface reflectance in wavelength interval, AA. 

Precise estimation of BBA from NBR from a remote 
sensing image can be obtained when following three factors 
are accounted for: (i) the spectral reflectance characteristics 
of the surface of interest; (ii) the spectral distribution of the 
irradiance; and (iii) the wavelength region of the discrete 
spectral bands. 

The AWiFS image has four narrow spectral bands in 
VIS (k = 520-590 nm and 620-680 nm), MR (770-860 nm) 
and SWIR (1550-1700 nm). The reflectance in these narrow 
bands (NBR) can be converted to the total BBA in the 
spectrum (400-2500 nm). Similarly, seven narrow spectral 
bands of MODIS image in VIS (459-479 nm, 545-565 nm 
and 620-670 nm), MR (841-876 nm and 1230-1250 nm) 
and SWIR (1628-1652 nm and 2105-2155 nm) can be used 
for NBR to BBA conversion. To convert from NBR to BBA, 
however, an improved understanding of the relationship 
between the reflectance characteristics of the narrow bands 
with that of broad band for snow surfaces is required. 

In the present paper, the spectrum was divided into 
different sub-regions: for AWiFS image the VIS sub region 
(400-700 nm); MR sub region (700-1500nm) and SWIR 
sub region (1500-2500 nm). In case of MODIS image, 
different sub regions are VIS (400-500 nm, 500-600 nm and 
600-700 nm), MR (700-1000 nm and 1000-1500 nm), and 
SWIR (1500-2000 nm and 2000-2500 nm). The selection of 
the sub regions depends primarily on: (i) consistency of data 
in AWiFS and MODIS narrow spectral bands in these 
regions; (ii) independence in spectral reflectance 
characteristics of snow surface in these regions. It is very 
important here to examine the relationships between AWiFS 
and MODIS spectral bands and their corresponding sub 
region as mentioned above. 

The total BBA for AWiFS image considering solar 
irradiance in the three sub-regions (i.e., VIS, MR, and 
SWIR) can be represented as: 
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a AWiFS 



TiAAl EsAAl P^xi + ^AA2 E sAA2 P^Xl + ^AA3 E sAA3 P^3 
Y*AA E sAA 



(6) 



where the wavelength interval AX represents spectral 
region (400-2500 nm); AX1 , A22 and AA3 correspond to 
three sub regions of the spectrum i.e. A/ll: 400-700 nm; 
AA2 : 700-1500 nm and AA3 : 1500-2500 nm. E s and p are 
solar irradiance and reflectance respectively in different 
bands. 



Similar to (6), the BBA for MODIS can be represented 



as: 



a M0DIS 



s?AXm 7 „ 

LAXml^sAXm PAXm 
T,AA E sAA 



(7) 



where AA ml _ m7 (m suffix for MODIS) represents seven sub 
regions of solar spectrum in VIS, NIR and SWIR. E sAAm an d 
p A0 1 „ are the solar irradiance and reflectance 

' HATH 1 — 771 / 

corresponding to broadband sub regions in different 
wavelength bands, AA ml _ m7 (AA ml > AA m2> AA m7 ) 
respectively. 

The estimated values of broadband reflectance obtained 
from (6) and (7) may be different from the narrow band 
reflectances in corresponding AWiFS and MODIS spectral 
bands respectively. Therefore, an adjustment factor that 
defines the relation between broad band and narrow band 
reflectance may be derived [30] . The adjustment factor 
quantifies the ratio of the reflected radiations within the 
narrow spectral bands to the reflected radiations of the 
corresponding broad band in VIS, NIR and SWIR regions. 
Let fi to f 3 and f lm to f 7m be the adjustment factors for 
AWiFS and MODIS images respectively. Equation (6) for 
AWiFS can then be rewritten as: 



_ fl 1 E sAAvPaAv +h 1 E sAAn P^An + h 1 E sAAsP^As 



T.aa e s, 



(8) 



where subscripts v, n and s refer to narrow VIS, NIR 
and SWIR bands of AWiFS sensor respectively. By 
comparing (6) with (8), f x , f 2 and f 3 can be estimated for 
AWiFS as: 



/' Pav E sAAv 

l — ^7 

L t sAAl PaAI 



h = 



PAn E sAAn 
Z E sAA2PaA2 
PAs E sAAs 



T, E sAMPaA3 

Similarly, the adjustment factors, f lm , f 2ri 
may be estimated for MODIS as: 

PAvl E sAlvl 



flm 
flm 
hm '■ 
Urn ] 
fsm 
fbm 



E E sAAm 1 PAXm 1 
. PAnl E sAXnl 

Y* E sAXm2PAXml 
PAv3 E sAXv3 

T. E sAAm3 PAXm3 
PAv4 E sAAv4 

X E sAhn 4 PAkm 4 
, PAnS E sAlnS 

?. E sAZm5 PAAmS 

PAs6 E sAAs6 
£ E sAXm 6 PAXm 6 



(9a) 
(9b) 
(9c) 

> flm > 

(10a) 
(10b) 
(10c) 
(lOd) 
(10e) 
(100 



fi 



7m 



PAs7 E sAAs7 



(10g) 



X E sAlm 7 PAXm 7 

where subscripts v, n and s refer to narrow VIS, NIR and 
SWIR bands of MODIS sensor. AX ml , AX m2 AX ml 
represent broadband sub-regions consistent with MODIS 
spectral bands. 

Assuming a linear relationship between the narrowband 
and the broadband albedo [2, 9] , the surface albedo from data 
of AWiFS and MODIS sensors can be estimated from: 



a AWiFS = C 1 p 2 + C 2 p 4 + C 3 p s +d 
dMODIS =C lm p lm + C 2m P 2m + C 3m p 3m +c 4m p 4m 

+C5m Psm +C 6m P 6m +Clm P 7m +d m 



(ID 

(12) 



where c's are coefficients, p's are narrow band reflectance 
in AWiFS and MODIS bands, and d is an offset, which is 
typically determined by comparing estimated values of 
albedo from model with those obtained from in-situ 
observations. On combining (8), (9) and (11), the values of 
coefficients may be estimated as: 

c = u--T L (13) 
where '.i's are the percentages of solar radiation within 
the wavelength bands of AWiFS and MODIS channels 
considered respectively. This is estimated using in-situ 
observations of solar spectral irradiance (Fig. 6). The values 
of / are calculated using the in-situ observations of spectral 
reflectance and spectral irradiance in (9) and (10) recorded 
at the time of satellite pass. 
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Fig. 6 Field observations of spectral reflectance (0-1) and transmitted 
solar spectral irradiance for different dates of AWiFS acquisition 
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VII. RESULTS AND DISCUSSION 

A. Spectral Reflectance and Spectral Irradiance 

The upwelling 't v ' and downwelling V transmissivity 
due to Rayleigh scattering from ground to sensor and sun to 
ground respectively for AWiFS and MODIS bands are 
given in the Table IV. Fig. 6 shows the average spectral 
reflectance characteristics of snow surface and spectral 
distribution of transmitted solar irradiance at different 
sampling locations surrounding Solang and Dhundi during 
six different passes of satellite carrying AWiFS sensor. It 
can be inferred from these figures that snow is characterized 
by high spectral reflectance in the visible part of an 
electromagnetic spectrum. It tends to decline in the near 
infrared region until 1050 nm wavelength, where marginal 
increase in reflectance occurs. A minor peak in reflectance 
can be seen at approximately 1090 to 1100 nm wavelength. 
The minor peaks can also be seen at around 1800nm and 
225 Onm wavelength with a sharp dip in reflectance at about 
1950 nm to 2050nm. This is due to the presence of water 
absorption bands at these wavelengths. The transmitted 
solar irradiance (Figs. 6a-c) has been found to be the 
maximum in visible wavelength region with a gradual 
decrease with increase in wavelength. The irradiance 
contribute 46.2 % in 400-700 nm, 26.8% between 700-1000 
nm, 16.7% in 1000-1500 nm and 10.2% in 1500-2500 nm 
wavelength regions. 

These in situ observations have also been used to: 

i) establish a relationship between existing 
narrow AWiFS and MODIS wave bands and corresponding 
broadband for estimation of values of the conversion factor 
/using (9) and (10); 

ii) estimate BBA and its validation with AWiFS 
satellite sensor albedo. 



The values of empirical coefficient, C, for AWiFS and 
MODIS sensors have been estimated using (13) and are 
given in the Table V. 

TABLE V MODEL COEFFICIENTS FOR AWIFS AND MODIS TO CONVERT 
NARROW BAND TO BROAD BAND ALBEDO 



AWiFS Coefficient 


MODIS Coefficient 


Q 


0.463 


Clm 


0.145 


c 2 


0.360 


C2m 


0.275 


c 3 


0.094 


C3 m 


0.138 


d 


0.026 


C4.n1 


0.165 






Csm 


0.214 






Cgm 


0.06 






C7m 


0.06 






d m 


-0.011 



B. Broadband Albedo ( BBA ) and Validation 

In this section, results of satellite sensor estimated snow 
BBA and their validation with in-situ observations at Solang, 
Dhundi and Patsio sites have been presented. The non-snow 
area in satellite images is masked using normalized 
difference snow index (NDSI >= 0.4) and reflectance in 
NIR > 0.10 [31] . Broadband albedo validation details at 
different sites are given below: 

1) Validation at Solang and Dhundi Sites Using Awifs 
Image: 

Figs. 7 (a-f, g-1) show the topographically uncorrected 
and corrected AWiFS reflectance images of Beas basin 
respectively. The thematic maps of (surrounding Solang and 
Dhundi) radiometrically corrected BBA obtained using (11) 
for snow covered regions are shown in the Figs. 7 (m-r). 
The value of radiometrically corrected broadband albedo of 
snow in these maps varies from 0.2 to 0.9. 
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Fig. 7 (a-f) Topographically uncorrected AWiFS images of different dates of Beas Basin (surrounding Solang and Dhundi); (g-1) Topographically 
corrected images; and (m-r) Thematic maps of broadband albedo for snow covered region after masking non-snow area 
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Fig. 8 Comparison of the retrieved snow broadband albedo with the in-situ measurements at sites Solang and Dhundi 
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In majority of the area, higher values (0.7-0.9) can be 
seen for the data acquired during January and February 
months than those obtained for the month of March, where 
the albedo in most of the snow area ranges from 0.5-0.7. 
This decrease may be attributed to the melting of snow 
during this month. The results on validation of satellite 
estimated BBA of snow using AWiFS with in- situ 
observations (data obtained using optical spectroradiometer) 
are shown in the Fig. 8. The overall R between observed 
and estimated values exceeds 0.94. The corresponding 
RMSE value has been determined as 0.03. These values 
indicate that the estimated BBA retrievals for AWiFS image 
bear a close relationship but may not be considered as 
identical to the in situ measurements (Fig. 8). This may be 
attributed to the mismatch in spatial resolution of the image 
56 m with in situ measurements, which happen to be point 
data. 

2) Validation at Patsio Sites Using Awifs and MODIS 
Data: 



The Patsio site lies in the middle Himalaya where the 
temperatures are quite low as discussed earlier. The snow 
BBA measured from albedometer on AWS has been used 
for validation of BBA estimated from proposed algorithm 
applied at Patsio site. Topographically uncorrected and 
corrected AWiFS reflectance images of lower and middle 
Himalaya (all three stations Solang, Dhundi and Patsio sites) 
have been shown in the Figs. 9 (a-g, h-n) respectively. Figs. 
9(o-u) show the thematic maps of snow BBA, which have 
been validated with in- situ albedo measurements at Patsio 
site. Fig. 10 shows the relationships between (AWiFS) 
retrieved snow broadband albedo with the in- situ 
measurements at Patsio site. The overall R~ between 
measured and estimated values for AWiFS data exceeds 
0.94 and the RMSE has been obtained as 0.0274. The R 2 
value and RMSE show that the AWiFS estimated BBA is in 
close agreement with the in-situ measurements at the Patsio 
site. 
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Fig. 9 (Top Row) Uncorrected AWiFS images of different dates includes all three sites Solang, Dhundi and Patsio; (Middle Row) Atmospherically and 
topographically corrected images; and (Bottom Row) Thematic maps of broadband snow albedo after masking non snow area 
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Fig. 10 Comparison of the satellite (AWiFS) retrieved snow broadband albedo with the in-situ measurements at site Patsio 
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Fig.l 1 (First Row) MODIS images of different dates includes all three sites Solang, Dhundi and Patsio; (Second Row) Thematic maps of broadband snow 

albedo after masking non snow area 



The snow BBA estimated from MODIS images at Patsio 
site has been validated with in-situ observations and global 
albedo model proposed by [11]. Figs. 11 (a-i) show the 
atmospherically corrected MODIS images of different dates. 
Corresponding thematic maps of BBA for snow covered 
region are shown in the Figs. 11 (j-r). Fig. 12 (a) shows the 
correlation between MODIS estimated BBA of snow using 
proposed model and in-situ observed albedo at Patsio site 
while Fig. 12 (b) shows the correlation between MODIS 
estimated BBA of snow using model and in-situ 
observed albedo at the same site. Fig. 12 (c) shows the 
RMSE of the proposed model and Liang model A high 
correlation (ft ~ = 0.88 for proposed model, and R~ = 0.84 
for Liang model) and low RMSE (0.026 for proposed model 
and 0.031 for Liang model) have been observed between 
estimated albedo and in-situ measured albedo for both the 
models although proposed model performed slightly better 
for the validation site. 



The performance of the proposed model has also been 
tested rigorously for a large snow covered region of 
Himalaya shown in Figs. 13 (a-b) vis-a-vis global albedo 
model proposed by [11] for nine images of different dates. 
A high value of correlation (ft ~ better than 0.99) and a low 
value of RMSD (0.015) have been observed between the 
two models (Fig. 14) for snow covered region. The global 
albedo model developed by [11] has not been validated for 
snow and ice, although the model is validated for other land 
covers. These results show that the Liang model is also 
performing well for snow and ice land covers. 
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Fig. 12 Validation of the satellite (MODIS) retrieved snow broadband 
albedo with the in situ measurements at site Patsio (a) Proposed Model (b) 
Liang Model (c) RMSE 



Fig. 13 (a) MODIS image of large snow covered Himalayan region and (b) 

Snow covered region after masking non snow area used for comparative 
analysis between proposed broadband snow albedo model and Liang global 

model 
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Fig. 14 Comparative and correlation analysis between proposed snow 
broadband albedo and Liang global model using different dates MODIS 

data 
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VIII. CONCLUSIONS 

Snow surface broadband albedo is a fundamental 
component in modeling any surface processes. In this paper, 
an operational algorithm for direct retrieval of snow 
broadband albedo using AWiFS and MODIS data has been 
developed and presented. This appears to be the first 
algorithm reported to estimate braodband albedo of snow 
from narrow band reflectances using AWiFs images in the 
Himalayan region. The in-situ measurements of spectral 
reflectance and transmitted solar irradiance have been used 
for estimation of conversion factor and coefficients required 
for development of snow broadband albedo using AWiFS 
and MODIS data. The results of BBA of snow retrieved 
from AWiFS and MODIS have also been validated with in- 
situ measurements at Solang, Dhundi and Patsio sites. The 
overall R 2 and RMSE values between in-situ observed and 
estimated BBA for AWiFS were observed better than 0.94 
and 0.03, while for MODIS data it were 0.88 and 0.026 
respectively. For MODIS data, the results from the proposed 
algorithm have also been found to be in agreement (overall 
R 2 , 0.99 and RMSD, 0.015) with published global albedo 
model for the large Himalayan snow cover region. 
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